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I.  Problem  Description 

Problems  of  so-called  decentralized  observation,  coordination  and 
hypothesis  detection  are  of  increasing  current  interest,  and  rather  little 
concrete  information  is  available  on  them.  Here,  we  discuss  a  simulation 
study  of  one  of  the  simplest  of  such  problems,  in  order  to  provide  some 
insight  into  the  phenomena  that  might  be  expected  in  more  general  cases, 
and  to  provide  a  guide  to  analysis.  We  present  these  results,  limited 
though  they  are,  because  of  the  importance  of  the  class  of  problems,  and 
the  paucity  of  available  information  and  intuition. 

We  deal  with  a  decentralized  detection  problem  involving  two 
observers,  denoted  X  ,  Y,  and  a  coordinator.  There  are  two  possible 
hypotheses  Hq  and  ,  with-  ^  *  P  {H  =  H^}  given.  Each  observer 

takes  an  observation  at  time  1  and  may,  if  it  wishes,  take  an  observation 
at  time  2.  Let  ^(X),  w.(Y)  denote  the  conditional  probability  of 

X,  Y,  respectively,  given  i  observations.  The  observers  do  not 
communicate  with  each  other,  and  at  the  end  of  their' observation  periods', 
each  transmits  its  w.  to  the  coodinator  who  then  computes  the  system 

posterior  probability  and  decides  on  either  or 

If  HQ  is  chosen  when  Hj  is  true,  there  is  a  cost  a0l  •  ^ 

Hj  is  chosen  when  HQ  is  true,  there  is  r  a1Q  .  We  view  Hj 

as  being  more  important  than  HQ  ,  and  set  aci  '  a10  •  If  at  least 
one  observer  takes  a  second  observation  and  Hj  is  true,  then  there  is  an 
added  delay  cost  .  Then  the  problem  is  just  a  decentralized  version 
of  the  standard  'disorder'  problem.  The  cost  of  delay  under  Hj  forces 
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each  observer  to  somehow  take  the  other  into  account  (in  deciding  whether 
to  take  a  second  observation).  There  is  no  second  guessing  in  the  setup, 
since  there  is  no  implicit  communication  among  X,  Y  .  This  is  deliberate, 
since  even  the  behavior  of  the  described  structure  is  not  understood.  Define 
g(7r)  =  min[ajjjir  ,  a^Cl-ir)]  .  Then  the  sample  cost  is 

Cost(Tr)  =  g(n)  +  CqU  I  {at  least  one  observer  takes  a  second  observation}. 

The  object  here  is  to  gain  some  insight  into  the  way  in  which  each 
observer  accounts  for  the  existence  of  the  other;  e.g.,  how  the  decision 
regions  differ  from  those  in  the  one  person  case,  the  parametric  dependence 
and  sensitivity  of  the  regions,  to  check  conjectures  concerning  qualitative 
behavior,  bounds  provided  by  one-person  cases,  and  to  check  for  'unusual' 
behavior  such  as  the  asymmetries  in  the  decision  regions  such  as  exhibited 
in  [1]  even  when  the  observers  had  the  same  distributions. 

II.  Background  Calculations 

Sufficient  statistics.  Let  Xj,  X2,  Y^,  denote  the  four 

possible  observations,  and  X1,  Y1  the  first  i  of  X,  Y,  respectively. 

We  suppose  that  {Xj,*  X^,  Y^,  Y^}  are  independent,  under  Hq  and  . 

Then 

\(X)  =  PCX1  |  H^PCHp/PCX1)  . 

If  X  takes  i  observations  and  Y  takes  j  observations,  then 
CD  *  -  P(Hj |  x\Yj)  =  P(X1,Yj  |  HjJPCHjJ/PCyV)  = 


a 
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=  PCX1  I  H1)P(Yj  I  HpPCHjJ/PCX1^) 

P(Hj  |Xi)P(H1  |Y^)P(Xi)P(Y^)/  ttq 
=  - 

l  PCH-IX^PCH L  |YJ)P(Xl)P(YJ)/P(H  ) 
k=0  K  k 

^(XJ^CY)/  % 

=  7Ti(X)YY)/TI0  +  (l-TT.(X)Hl-1T;j(V)y/(l-TT0) 


Thus  it  is  enough  for  X,  Y  to  transmit  their  conditional  proba¬ 
bilities  at  the  end  of  their  observation  period. 

The  choice  of  P(Xj|H^),  P(Y^|H^)  .  In  the  simulations,  these  did 

not  depend  on  j  ,  were  identical  for  X  and  Y  ,  and  we  choose  an  expon¬ 
ential  form.  Let  XQ  >  0,  X^  >  0.  Then  (density)  p(X^|H^)  =  X^exp  -  X^X^  . 

Since  up  to  4  observations  can  be  taken,  the  evaluation  of  the  cost  functions 
can  require  complicated  numerical  integrations  in  general.  With  the  expon¬ 
ential  distribution,  algebraic  expressions  can  be  obtained  for  these 
integrals,  and  the  computations  are  considerably  simplified.  Some  particular 
aspects  of  the  results  do  seem  to  depend  on  the  exponential  assumption,  but 
most  of  our  observations  and  remarks  probably  hold  much  more  generally. 


property  of  the  exponential  distribution.  We  have 


HpXjfexp  -  XjXp 

VX)  =  UqX^xp" -  XjXp  ♦  (l-rr0U0exp  -  XQX1) 
=  V i^ tv i  +  (1~Vxoexp  ~  *1 
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Let  •  Then 

(3)  yx) e  i°*  Vi/[Vi +  (l_W  5  active  region 

If  xo  >  \  ■ 

(4)  YX)  e  iVl/lVl  +  (1_VV*  ^  H  active  region 

The  active  region  depends  on  ^  and  the  X^  ,  and  is  never  the  entire 

segment  [0,1].  In  general,  some  ranges  for  the  posterior  probability -will 
be  relatively  unlikely  under  one  hypothesis  but  not  under  another,  so  we 
do  not  believe  that  the  exponential  distribution  is  too  specialized. 

Sufficient  statistics  for  Decision  (the  exponential  assumption  is 
not  needed  here.)  Based  on  X^ (respectively  Y^)  ,X (respectively  Y) 

decides  whether  to  take  another  observation.  We  show  that  it  is  sufficient  for  X 
to  use  ir^(X)  rather  than  the  observation  (and  similarly  for  Y's  decision). 

Let  d^,  dy  denote  the  decision  rules;  i.e.,  these  are  sets 

Y  Dy  such  that  (e.g.)  X  stops  if  X^  €  .  We  want  (with  hopefully  obvious 

notation) 

min  E_  cost(DY,Dy) 

Dx  ,  Dy  0  A 

(5)  *  min  E_  min  E_  Y  cost  (DV,D  ) 

Dy  0  Dx  *0,1  Y 

*  min  E_  min  < E.  v  cost  (X  goes,Dv), 

*0  <  %,X1  Y 


E  Y  cost  (X  stops, Dv) 

VA1  Y 


Note  that  (again  with  hopefully  obvious  notation,  and  using  (1)  to  write  g(n) 


in  terms  of  the  relevant  it.) 

i 


Ew  x  cost  (X  stops,  Dy)  = 

0,  1 


E^o  Xl  *CWL(X)  ,»1CY))ICY  €  Dy)  ♦  (c0  +  g(TT1(X),Tr2(Y))I(Y  e 


?7T0  ,Xj  (X2*Y1  ’Y2)  Ptt0  (X)  (X2  ,Y1  ,Y2) 


Thus,  (6)  is  a  function  of  no  and  itj(X)  only.  Similarly  for  the 

En  x  cost  (X  goes,  Dy)  .  Thus,  Dx  and  Dy  depend  on  ^(X),  ^(Y) 

0,  1 


not  on  Xj,  Yj  directly,  and  we  write  the  decision  regions  as  ^  and 


We  note  that  the  same  assertion  holds  in  the  N-step  problem. 


where  each  player  can  stop  and  communicate  to  the  coordinator  at  any  time 
T  N,  and  the  cost  for  continuing  to  observe  is  c^[time  that  the  last 

player  communicates].  The  proof  uses  a  simple'dynamic  programming  like' 
backward  induction,  and  is  omitted. 


The  calculation  of  the  decision  regions  .  The  above  cal¬ 

culations  imply  that  the  optimal  decision  regions  can  be  calculated  by 
fixing  ^yy  .  In  the  2-step  exponential  distribution  case,  the  optimum 

role  takes  the  following  form:  there  are  Tx  <  Tx  ,  Ty  <  Ty  such  that 
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such  that  the  second  X-  observation  is  taken  if  £  ^(X)  £  , 

and  similarly  for  Y  .  This  has  not  yet  been  proved  in  general,  but 
we  specialize  the  rules  to  this  class. 

In  our  simulations,  the  observations  of  X,  Y  had  the  same  dis¬ 
tribution,  under  either  or  .  By  the  comments  in  the  last  sub¬ 

section,  there  is  an  f(*)  such  that  given  Ty  =  (Ty,  Ty) ,  the  corres¬ 
ponding  optimum  satisfies  =  f(Ty)  .  By  symmetry,  Ty  =  f(T^). 

Thus  (*  denotes  optimum) 

(8)  T*  =  f(f(T*))  . 


In  the  calculations  (and  in  Figures  1,  2),  we  simply  calculated 
T^  given  Ty  by  checking  the  sign  of 


(9)  Ejj  ^(x)  (cost  if  x  stops,  given  Ty) 


5Tt  ,jt  (X)  (cost  if  x  8oes>  given  Ty)  . 
0*  1 


III.  The  Threshold  Tables. 

Before  presenting  and  discussing  the  numerical  results,  we  make  two 
remarks.  First,  the  thresholds  depend  on  the  lt^(X),1Ij(Y)  are  not 

sufficient  statistics  for  decision  in  the  usual  sense.  This  is  obvious 
from  (1).  But  the  ^-dependence  is  often  very  small  (see  Table  1).  Such 

a  dependence  does  not  occur  in  the  analogous  stan.dard  one  person  case. 

Even  with  ^  (X)  known,  still  provides  additional  information  on  the 
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distribution  of  the  ir  and  of  I{Y  goes}  .. 

Secondly,  the  lower  thresholds  were  quite  small  in  our  examples,  and 
the  upper  thresholds  (say  T^)  were  not  particularly  sensitive  to  variations 

in  Ty  between  0  and  the  optimum  .  So,  in  most  simulations,  we  simply 

fixed  Ty  at  a  small  value  -  in  that  range,  calculated  f(T^)  =  (T^.T^), 

then  optimized  over  Ty  .  Refer  to  Tables  1  to  4.  We  call  (XQ  =  a,A^  =  b) 
the  (a»b)  case. 


TABLE  1. 

Optimum  Upper  Thresholds  =  Ty  , • 

AQ  =  1,  Aj  *  3; (1,3)  case. 

(*) active  region  =  [0,.25],  always  continue 


TABLE  2. 

(3,1)  Case 

(**),  active  region  =  [.5,1],  always  stop. 


.25 


.5 


.73 


(1,2)  Case 


TABLE  4. 
(2,1)  Case 


Comments  on  the  tables.  In  all  crises,  the  optimum  thresholds 
for  X,  Y  were  the  same.  There  were  no  asymmetries  as  in  [  1] ;  this  is, 
in  fact,  one  of  the  phenomena  which  we  were  looking  for.  The  Cp- 
dependence  indicated  in  the  tables  is  as  expected.  The  ’C.-  dependence 

is  obvious,  but  for  most  of  the  range,  the  dependence  is  slight. 

The  Hq-  dependence  can  be  partly  explained  by  noting  that  (say,  from 
X's  point  of  view)  provides  information  on  ^(Y)  and  on  the  ul¬ 
timate  w  .  Loosely  speaking,  a  smaller  implies  that  w  will  be 

smaller,  for  fixed  ^(X).  This  ^p-effect  does  not  occur  in  the 


standard  one  person  problem. 

Owing  to  the  fact  that  the  'active*  region  is  a  proper  subset 
of  [0,1],  it  is  possible  (see  Tables  1,2)  that  the  optimum  decision  is 
to  always  stop  (or  to  continue)  for  all  attainable  values  of  the  ^(X) 

(or  ^(Y)),  for  the  given  ^  .  (This  'active  region'  effect  might  be 

exaggerated  for  our  case,  but  in  general,  it  simply  reflects  the  fact 
that  certain  7T^(X)-s.ets  are  relatively  unlikely  to  occur,  given  a  ^g.) 

In  the  analogous  one-person  problem  the  threshold  for  the  (1,3) 
case  is  greater  than  the  threshold  for  the  (3,1)  case*  There  is 'less 
risk’ in  continuing  in  the  former  case,  since  there  is  a  greater  smoothing 
of  the  decision  error  cost.  See  Figs. 3-5  and  further  comments  below. 

(Of  course,  in  the  one-person  case  the  thresholds  are  ^  -thresholds.) 

In  the  two-person  case  the  situation  is  reversed, in  the  cases  studied.  See 
the  tables.  We  have  not  yet  found  a  completely  satisfactory  explanation, 
but  it  might  be  partially  explained  by  the  fact  that  the  same  value  of  tt^(X) 

in  the  (3,1)  and  (1,3)  cases  implies  a  smaller  X^  in  the  (3,1)  case. 

Similarly,  the  same  observation  implies  a  larger  ^  (X)  in  the  (3,1) 

case.  The  reversal  cited  above  would  not  occur  if  the  thresholds  for  X 
were  written  in  terms  of  X^  (and  similarly  for  Y) .  These  remarks  point 
out  that  the  use  of  one-person  formulations  in  order  to  gain  insight  for 
the  two-person  case  must  be  done  with  great  care.  Similar,  though  not 
so  extreme,  observations  hold  in  the  (2,1),  (1,2)  cases.  See  Tables  3,4. 

An  additional  reason  for  the  'reversal'  is  the  fact  that  the  active  regions 
in  the  two  cases  start  at  opposite  ends  of  the  [0,1]  interval. 


Let  the  one^person  case  thresholds  be  (TL,  T^)  .  Then 

for  >  A^,  the  simulations  yielded  T^  >,  T^  =  T^  ,  T  _>  “  Ty  . 

The  reverse  occured  when  A^  >  A^  but  the  observation  suggests  that 

useful  information  might  be  obtainable  from  one-person  case  results  in 
special  cases. 

IV.  The  Threshold  Curves  -  Figures  1,2. 

By  the  argument  associated  with  (5) ,  the  optimum  thresholds  can 
be  obtained  by  fixing  T^,  Ty  ,  optimizing  over  T^,  T^  ,  and  then  over  the 

first  pair  by  seeking  the  fixed  points  of  (8).  As  noted  above,  the  lower 
thresholds  were  small,  and  varying  them  in  the  interval  near  zero  in  which 
they  normally  ranged  made  little  difference  on  the  obtained  values  of  the 
upper  thresholds  -  so  to  get  the  plotted  curves  we  set  T^  to  some  value 

U 

within  its  normal  range  and  varied  only  T^  .  Except  for  Fig.  la,  only 

the  corresponding  optimum  T^  are  plotted.  The  T^  were  calculated  by 
checking  the  sign  of  the  difference  between  the  arguments  of  the  inner  minimum 
on  the  right  side  of  (5) .  The  flattening  of  the  right  part  of  the  curve 
in  Fig.  1  is  due  to  the  value  of  the  active  region.  In  Fig.  la,  for 
example,  the  threshold  value  1/2  could  be  replaced  by  any  value  1/2, 
since  [0,1/2]  is  the  active  region.  In  Fig.l,  the' important'  hypothesis 
Hj  is  the  more  difficult  to  detect. 

By  (8),  only  the  points  (1),  (2),  (3),  in  Fig.  la  could  be  optimal 
pairs.  In  fact,  only  (1),(3)  can  be  optimum,  and  (2)  is  a  saddle.  Point 
(1)  is  optimum.  In  all  cases  tested,  the  optimal  pairs  were  on  the  diagonal. 
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This  will  always  be  the  case  if  never  decreases  as  Ty  increases. 

Since  the  ’continuation’  cost  must  be  paid  if  either  observer  continues, 
the  ’monotonic'  property  probably  holds  fairly  generally.  In  Fig.  1  note 
the  relative  insensitivity  of  the  optimal  to  Ty  (especially  as 

compared  with  the  Fig.  2  case).  This  is  connected  with  the  lower  threshold 
in  the  first  case.  As  increases  this  sensitivity  increases,  and  it 

decreases  as  ^  increases  (loosely  speaking,  as  ^  increases,  X 

’increasingly  ignores’  Y) .  In  Fig.  1,  compare  the  optimal  T^  with 

the  optimal  when  Y  does  not  take  a  second  observation  (this  is 

the  vertical  intercept).  It  is  virtually  equal  to  the  optimal  threshold; 
each  observer  essentially  ignores  the  others  possible  second  observation. 
Thus  useful  simplifying  rules  might  be  found  in  special  cases. 

Now,  refer  to  Fig.  2.  In  Figs.  2a, 2b,  the  optimal  thresholds  are 
unique.  The  optimum  T^  is  relatively  insensitive  to  modest  variations 

in  Ty  .  The  sensitivity  appears  to  be  greater  in  Fig.  2c,  which  also 

illustrates  how  complicated  the  situation  can  get.  However,  the  obser¬ 
vations  in  the  next  section  below  suggest  that  even  here  the  cost  is 
relatively  insensitive  to  modest  variations  in  T^,  Ty  about  their  optimal, 

or  to  variations  within  the  interval  where  the  curves  in  Fig.  2c  stay  near 
the  diagonal. 
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V.  The  Decision  Cost  Curves  -  Figure  3. 

Legend  for  Fig.  3.  Curve  (0)  denotes  gO)  =  min  [aQ17r  .a^Cl-w)  ] , 
curves  (i) ,  i  =  1,2,3,  indicate  the  average  decision  costs  for  the  one-person 
case  where  either  1,2,  or  3  observations  can  be  taken.  Then  the  abscissa  is 
Curve  (S)  denotes  the  decision  cost  part  (not  the  part)  of  the 

second  argument  of  the  inner  min.  on  the  right  side  of  (5)  (under  stop) ,  and 
Curve  (G)  denotes  the  decision  cost  part  of  the  first  argument  (under  con¬ 
tinue)  .  l.e.,  the  curves  denote  the  average  decision  costs  for  both  decisions 

from  X's  point  of  view.  Ty,  Ty  (i.e.,  Dy)  were  fixed  near  their  optimum 
values.  For  these  S,G  curves  the  abscissa  is  Ttj(X),  and  is  fixed. 

In  all  cases  Curve  (2)  <_  Curve  (S)  <_  Curve  (1) ,  Curve  (3)  <_  Curve  (G) 

£  Curve  (2).  This  is  intuitively  reasonable,  and  can  be  proved  to  hold  in 
general.  Since  Curve  (S)  is  the  decision  cost  conditioned  on  (Xj,tTq),  it  is 
averaged  over  Y's  one  or  two  observations.  Similarly  Curve  (G)  is  the  cost, 
averaged  over  Y's  one  or  two  observations  and  an  additional  X-observation. 

As  Up  increases  Curve  (G)  +  Curve  (2)  and  Curve  (S)  +  Curve  (1) 
(also  a  general  result),  since  as  ^  increases  the  information  (the 

additional  information  that  Y  provides  the  coordinator,  from  X's  point 
of  view)  in  Y's  possible  future  observations  decreases.  These  observations 
might  provide  some  help  in  simplifying  solving  for  the  thresholds  in  more 
complicated  cases. 

VI.  Total  Cost  Curves  -  Figs.  4,5. 

Figures  4,5  plot  expression  9  with  Ty,Ty  fixed  near  their  optimum 


values. 


Legend.  Curve  (T)  =  our  2  person  case,  with  abscissa  =  ^(X),  and  Ty/Iy 
fixed  near  their  optimum  and  ir^  fixed.  Far  i  =  1,2,  Curve  (i)  =  average 
cost  for  no  observations  minus  average  cost  given  i  additional  observations. 
Here  w^(X)  is  irrelevant  and  the  abscissa  is  Curve  (S)  is  as  Curve  (T) , 

but  with  Ty  =  Ty  a  0  (Y  never  continues).  Curve  (R)  is  as  Curve  (S) ,  but 
X  takes  2  additional  observations  if  it  elects  to  continue.  The  curves  are- 
typical  of  the  runs  taken.  The  zero  crossings  are  the  optimal  X-thresholds, 
for  the  given  Y-thresholds.  In  Fig.  4,  the  T,S,R,1  curves  have  similar 
structure,  the  sensitivity  at  optimum  being  slightly  less  for  the  T  curve. 

Now  refer  to  Fig.  5  for  the  (3,1)  case.  Here,  the  T-curve  is 
much  less  sensitive  to  ^(X).  The  nature  of  the  curve  implies  that 

modest  variations  in  the  thresholds  will  affect  the  total  average  cost 
only  slightly.  This  appears  to  be  the  case,  for  situations  such  as  in 
Fig.  2c,  when  there  are  multiple  crossings  of  the  diagonal.  In  such 
cases,  it  appeared  that  the  total  average  systems  cost  did  not  vary  too 
much  as  the  thresholds  varied  on  the  interval  between  the  local  minima. 

t 

Of  course,  these  remarks  are  merely  suggested  by  curves  such  as  those  in 
Fig.  5.  and  are  not  mathematically  proved  assertions.  Nevertheless,  they 
do  suggest  that  calculating  good  thresholds  might  not  be  too  hard  even  in 
cases  which  lead  to  curves  such  as  that  of  Fig.  2c. 
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FIGURE  4  TOTAL  COST  CURVES. 
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FIGURE  5  TOTAL  COST  CURVES. 


